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Abstract. We performed numerical simulations of stel- 
lar occultations by extra-solar cometary tails. We find 
that extra-solar comets can be detected by the appar- 
ent photometric variations of the central stars. In most 
cases, the light curve shows a very peculiar "rounded tri- 
angular" shape. However, in some other cases, the curve 
can mimic a planetary occultation. Photometric variations 
due to comet occultations are mainly achromatic. Never- 
theless, if comets with small periastrons have smaller par- 
ticles, these occultations could be chromatic with a larger 
extinction in the blue by few percents. 

We also estimate the number of detections expected in 
a large photometric survey at high accuracy. By the ob- 
servation of several tens of thousand of stars, it should be 
possible to detect several hundreds of occultation per year. 
We thus conclude that a spatial photometric survey would 
detect a large number of extra-solar comets. This would 
allow to explore the time evolution of cometary activity, 
and consequently would probe structure and evolution of 
extra-solar planetary systems. 

Key words: Occultations - Comets - Planetary systems 



1. Introduction 

The search for extra-solar planets through photometric 
variation is a well-known problem analyzed in many as- 
pects for several years (Schneider et al. 1990; Schnei- 
der 1996). However, planets are not the only objects 
detectable by their effect on stars' brightness, comets 
can also induce photometric variations. In our study of 
(3 Pictoris photometric variations (Lecavelier des Etangs 
et al. 1995), we concluded that the variations observed on 
November 1981 could be due to the passage in front of 
the star of either a planet or a dusty cloud (Lecavelier des 
Etangs et al. 1997, Lamers et al. 1997). If the latter is the 
correct explanation, this cloud of dust must be a cometary 
tail for two reasons. First, because the lifetime of the dust 
in such a system is very short and one must find a way to 
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produce the dust continuously. Second, because the cloud 
shape is constrained by the observed light curve; the cloud 
cannot be spherical, it should have a sharp edge in the part 
pointing toward the star and a huge cloud of dust in the 
opposite direction: exactly like a cometary tail (Lecavelier 
des Etangs, 1996). 

In the case of the solar system, the observations of stel- 
lar occultation by comets have been discussed by Combes 
et al. (1983). Observations of extinction and polarization 
of star light by dust of cometary tails in the solar system 
have been published (Dossin 1962, Ninkov 1994, Rosen- 
bush et al. 1994). 

But the solar system is only one particular planetary 
system at a given age. It is known that cometary activity 
was formerly much more important, and the well-known 
case of (3 Pictoris shows that during the first 10 s years, a 
planetary system is expected to show large cometary ac- 
tivity (Ferlet & Vidal-Madjar 1995, Vidal-Madjar et al. 
1998). Moreover, presence of comets around stars can be 
considered as an indirect signature for the presence of 
gravitational perturbations, and possibly caused by plan- 
ets. 

In this paper, we deal with the important possibility 
of detecting cometary activity from a photometric sur- 
vey. We take the COROT space mission as an example 
of what will be achievable in the very near future (Baglin 
et al. 1997). COROT, which primary aim is stellar seis- 
mology, will be launched in early 2002. It will allow a 
survey of about 30 000 stars with a photometric accuracy 
of few 10~ 4 during several months. Here, we predict the 
probabilities of detecting comets using such a photometric 
survey. 

In Sect |^, we describe the model of the cometary occul- 
tation, and then give the expected light curve in Sect | 
Estimates of the number of comets which could be de- 
tected are given in Sect |[ The conclusion is in Sect |[ 

2. A model of cometary occultation 

As pointed out by Lamers et al.(1997), a geometric dis- 
tribution of the dust must be assumed in order to eval- 
uate the photometric stellar variation due to extra-solar 
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comets. Then, by taking into account the optical proper- 
ties of the cometary grains, the photometric variation can 
be estimated. 



2.1. Distribution of dust in a cometary tail 

The evaluation of the distribution of the dust in a 
cometary tail can be made through particle simulation. 
The input parameters are the comet orbit, the dust pro- 
duction rate, the ejection velocity and the size distribu- 
tion. We assume a size distribution dn(s) of the form 



dn(s) = 



(i - s / s y 



(i) 



as observed in the solar system, where s is the dust size. 
We take sq = 0.1pm, n = 4.2, m = n(s p — so)/sq, and 
s p — 0.5 pm (Hanner 1983, Newburn & Spinrad 1985). 
This distribution starts at so = 0.1/im with dn(so) = 0, 
peaks at s p — 0.5/iin, and has a tail similar to a s~ n 
distribution for large sizes. 

Dust sensitivity to radiation pressure is given by the ft 
ratio of the radiation force to the gravitational force. We 
take 



(3 = 0.2 



»/M. \ 



Lq/MqJ \lpm 



(2) 



where L*/M* is the luminosity-mass ratio of the star. So 
the orbits of the small grains are more affected by radi- 
ation pressure than those of large grains. This is a very 
good approximation for particles larger than 0.1/im of any 
realistic composition and for the solar spectrum (Burns et 
al. 1979). 

The dust production rate P is assumed to vary with 
r, the distance to the star, and is taken to be 



P = Po 



(3) 



(see, for example, A'Hearn et al., 1995; Weaver et al., 1997, 
Schleicher 1998). The dust production is taken to be zero 
beyond 3 AU for a solar luminosity. 

As soon as a grain is produced from a comet nucleus, 
we assume that it is ejected from the parent body with 
a velocity u e ject in an arbitrary direction. The grain then 
follows a trajectory defined by gravitation and radiation 
forces. The ejection velocity depends upon the particle 
size. We take v C jcct = {A + Byf/3) (Sekanina & Lar- 
son 1984), which approximates the results of Probstein's 
(1969) two-phase dusty-gas dynamics for the acceleration 
by the expanding gas within tens of kilometers from the 
nucleus. The coefficient A and B depend upon many pa- 
rameters such as the thermal velocity of the expanding 
gas. We used A = B = 1 s km -1 which is a good ap- 
proximation of different values measured for the comets 
of the solar system (Sekanina & Larson 1984; Sekanina 
1998). This gives an ejection velocity of rs 585km s _1 



for s = 0.1/im. The ejection velocity is smaller for larger 
grains which have smaller cross section area to mass ratio 
(deject (10/im) « 124km s _1 ). We checked that any other 
realistic values for A and B give similar light curves within 
few percents. 

2.2. Stellar parameter 

The simulations have been performed with the mass, lu- 
minosity and radius of the central star set to solar values 
(A/* , and i?» ) . Simulation with other parameters could 
be possible. A larger mass for the central star would induce 
a shorter time scale; a larger luminosity would increase the 
effect of the radiation pressure on grains; a larger radius 
would decrease the relative extinction (Eq. ^). However, 
the conclusions do not strongly depend on these stellar 
properties. For instance, adopting the properties of an 
A5V star does not change the shape of the light curves. 
By scaling the production rate to the star luminosity, we 
found a quantitative change by less than a factor of two. 

2.3. Grain properties 
2.3.1. Extinction 



is 

ext IS 



The extinction cross-section of a dust grain 
Qext( s > A)7rs 2 , where the extinction efficiency, Q 
slightly dependent on the particle size (s) and radiation 
wave length (A). If <5 sca is the scattering efficiency (see 
Sect. |2.3.2 ) and Q a bs is absorption efficiency, Q ex t{s, A) = 
Qsca + Qabs- We take Q abs = 1 if s > A, and Q abs — s/X 
if s < A, which is a good approximation for optical wave- 
lengths and grains larger than 0.1/im (Draine & Lee 1984). 

The total extinction is hence calculated by adding the 
extinction due to all particles in the line of sight to the 
star. The optical depth r due to the dust is 



Spart. -^grain/part. Qext(s, A)7TS 2 



(4) 



where S is the projected area of the line of sight. In the 
simulation, one particle represents an optically thin cloud 
of several dust grains. The number of physical grains per 
particle in the simulation, N gIa x n / part , is 



N, 



grain/part. 



3M dust J dn(s) 
4:npN pav t. J s 3 dn(s) 



(5) 



where M^ust is the total dust mass, p is the dust density, 
and A^part. is the number of particles in the simulation. 
-Afpart. is set to a few 10 4 in order to keep reasonable com- 
puting time. 

Because the cometary cloud is optically thick but its 
size is smaller than the size of the star, we mapped the 
stellar surface through a set of cells in polar coordinates. 
For each cell i, we calculate the optical depth due to the 
particles within this cell of area Si = wR*)- The 
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ratio of the flux observed through the cloud (F ext ) to the 
initial stellar flux (F*) is 



Fext v - ^ Si e 



TlRl 



(G) 



The number of cells is the best compromise between 
the spatial resolution and the number of particles in the 
simulation. We take into account that the maximum con- 
tribution to Tj by each particle must be < 10 _1 in order 
to achieve an accurate result in spite of the quantization 
of the extinction. Si and iVp ar t are thus constrained by 



SiN part 3M dust Q ext J s 2 dn(s) 



10 



Ap J s 3 dn(s) 



(7) 



The limb darkening is not taken into account in this 
work, because its effect is negligible. 

2.3.2. Scattering 

As already pointed out by Lamers et al. (1997), the main 
part of an occulting dust cloud is seen through a very 
small scattering angle. Thus, the total star light forward- 
scattered to the observer can be large because the phase 
function is strongly peaked to small angles for which the 
diffraction has the dominant contribution. The phase func- 
tion for the diffraction is [] 

2 



T.2 



P(6,x) 



47T 



Ji (x sin ( 



x sm t 



(8) 



where x = 2tts/ X. The phase function has been normalized 
by J P(6)dn — 1 and depends upon the wavelength A. 
Examples of such phase functions are shown in Lamers et 
al. (1997). 

The scattered light (F sca ) is evaluated by adding the 
contribution of each particle in the simulation. 



= N, 



grain/part. / 

part. 



Q sca ns 2 P(6) _ 



(Ti„+T 



(9) 



Ti n is the total extinction along the path from the star to 
the scattering grain and r ou t is from the grain to the ob- 
server. The grain is at distance r from the star. For s > A, 
the scattering efficiency Q sca is assumed to be the diffrac- 
tion efficiency for large grains: Q SC a = Qd = 1 (Pollack & 
Cuzzi, 1980). For small particles (s < A), we used the basic 
approximation Q sca — (s/A) 4 (van de Hulst, 1957). The 
result is very insensitive to this last assumption because 
forward scattering at small angle is largely dominated by 
diffraction on particles larger than the wavelength. Be- 
cause each cloud represented by a particle is very thin, 
and the total extinction is small (rm, T out < 10 _1 ), we 
assumed single scattering (no source function). 



1 The Eq. 12 in Lamers et al. (1997) has a typo-mistake: x 
in the first fraction must be x 2 . 



In fact, for very peaked forward-scattering function on 
very small scattering angle, the finite size of the star must 
be taken into account (especially when the dust cloud is 
seen superimposed on the star surface) . Hence, we mapped 
the stellar surface by small arcs centered on the dust par- 
ticle. We then calculate the total scattering by adding the 
contribution of each arc. 

3. Light curves 

The light observed at a given time is the sum of two op- 
posite effects: the increase of the brightness through the 
scattering by particles at small angle from the line of sight 
and the decrease of the brightness through the extinction 
(including absorption) by particles on the line of sight. 

3.1. First order estimations 
3.1.1. Extinction 

From Eqs. || and ||, in the case s A (Q ex t — 2), and 
neglecting the inner extinction, we find an upper limit for 
the extinction: 



exp 



GMdust J s 2 dn(s) 
AirpR 2 J s 3 dn(s) 



(10) 



For 
(Sect 



the particle size distribution given by Eq. 



2.1), we have (/ s 2 dn(s))/(J s 3 dn(s)) = 0.15/z 



J 



Adopting a dust production of 10 6 kg s 1 during 100 days, 
we find M dust a 10 13 kg. Hence (AF/^) max a 10~ 3 . This 
order of magnitud e is c lose to the result of more elaborated 
calculations (Sect. 3^). We can already conclude that this 
variation could be detected by a space photometric survey. 

We note that the photometric variations of f3 Pic in 
1981 was ~ 5 • 10~ 2 . The above calculation gives an es- 
timate of the lower limit of the needed production rate : 
P k, 5 • 10 7 kg s -1 . This is consistent with P £ 10 8 kg s _1 
obtained from the estimation of Lamers et al. (1997): 
P/Av > 10 8 kg km -1 and with an assumed relative ve- 
locity between the dust and the nucleus of the comet 
Av 1km s^ 1 . 

3.1.2. Scattering 

If we take a comet at r = 1 AU from the star and an 
impact parameter on the line of sight of one half the stellar 
radius, then 9 w 7'. For the particle size distribution given 
by Eq. [j], the phase function is P{9 w 7') sa 700. Therefore, 



F sca t _ 3M dust P{6) J s 2 dn(s) 
F* Apr 2 J s 3 dn(s) 



4-10" 



(11) 



As a first conclusion, it is clear that both extinction 
and scattered light can be detected by very precise photo- 
metric survey, although scattering gives photometric vari- 
ations an order of magnitude smaller than extinction. Ex- 
tinction appears to be the major process observed when 
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the comet is passing in front of the star. This event can 
thus be called an occultation. 

3.2. Light curves 

Taking into account the dust spatial distribution and the 
extinction within the tail, we calculated the light curves of 
a set of cometary occultations. From a given comet orbit, 
and dust production rate, we compute the full motion of 
each dust particle; the result is the variation of the star 
light as a function of time. 

Two typical light curves resulting from the simulation 
are shown in Fig. [l] and|^. The majority of cometary occul- 
tations gives light curves with a very particular "rounded 
triangular" shape (Fig. Q). This occurs when the dense 
cometary head first occults the star and gives a very fast 
and sharp brightness decrease. It is then followed by the 
tail which gives an additional slow decrease. Similarly, 
the subsequent brightness increase is also sharp when the 
cometary head is going out of the occulting part. Then 
the increase slows down and the brightness returns to the 
normal level when the cometary tail is less and less dense 
in front of the star. 

However, in some configurations, the tail can be 
aligned with the line of sight. In these cases the light curves 
are more symmetric (Fig. ||). They can mimic planetary 
occultations (Fig. ||). Because of the noise, it will be dif- 
ficult from such observations to differentiate between a 
comet and a planet. 

3.3. Color signatures 



As seen in Sect. 2.3, the light variations show some color 
signature due to the optical properties of the grains. Par- 
ticles with a size smaller than the wavelength are less ef- 
ficient for extinction. Thus, extinction is smaller at larger 
wavelengths. The forward scattering is more peaked to the 
small angles at shorter wavelengths (Lamers et al. 1997). 
The scattering is thus larger at smaller wavelengths when 
the cometary cloud is occulting the star (a contrario, the 
scattering is larger at larger wavelength when the comet 
cloud is seen far from the star; but then the scattered 
light becomes negligible). As a result, there is a balance 
between the additional scattered light and the extinction; 
both are larger at smaller wavelengths. It is difficult to 
predict what will be the color effect on cometary occulta- 
tion light curves. 

In addition, most of the extinction and scattering due 
to the cometary dust is concentrated in the inner coma, 
where the number of large particles is larger because 
smaller particles are more efficiently ejected by radiation 
pressure (Sect. pTl] ) . As a consequence, in the optical, and 



with the size distribution assumed in Sect 2.1, the light 
curves are barely dependent on the wavelength (or color 
band) . The difference between variations in blue and in red 
is smaller than few percents. This color signature would 
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Fig. 1. Plot of the photometric variations during a 
cometary occultation in red (A ~ 8000 A, thick solid line) 
and in blue (A ~ 4000 A, thick short-dashed line). The 
insert is a view from the top when the comet is crossing 
the line of sight (Y = 0) at the periastron. The produc- 
tion rate is 2 • 10 6 kg s _1 at 1 AU. The scattered light is 
given by the top thin lines, the extinction by the bottom 
thin lines (long-dashed line for A ~ 8000 A; dotted line for 
A ~ 4000 A) . The total variation is plotted with the thick 
lines. This light curve presents the very specific "rounded 
triangular" shape observed in the majority of simulations 
of cometary occultations. The difference between the vari- 
ation in the blue and in the red is negligible and less than 
1%. 



be very difficult to observe. This is beyond the today tech- 
nical feasibility. 

However, for comets at small distances from the Sun 
0.5 AU), its appears that the dust size distribution is 
peaked at smaller sizes (Newburn & Spinrad 1985). This 
decrease of a p with the heliocentric distance could be due 
to particle fragmentation. Although it is difficult to guess 
the size distribution for extra-solar comets, it is easy to un- 
derstand that if small particles are more numerous (with 
s < 0.5/im), the occultation will show a larger color sig- 
nature. We checked that a comet with a periastron at 
0.3 AU and s p = 0.25^m (Newburn & Spinrad 1985) gives 
larger extinction by about 15% in the blue than in the red 
(Fig. |). 

3.4. Conclusion 

Photometric variations due to the stars' occultation by 
extra-solar comets could be detectable by photometric 
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Fig. 2. Same as Fig.l. Here all the parameters are the 
same except that the longitude of the periastron is at 
90° from the line of sight. This gives a production rate of 
5 • 10 5 kg s _1 at 2 AU. The cometary tail is aligned with 
the line of sight; consequently the light curve presents a 
more symmetrical shape resembling a planetary occupa- 
tion. The net color effect is hardly detectable and of the 
order of 1.5%. 



measurements with an accuracy of ~ 10~ 3 - 10 -4 . In 
most cases, the particular "rounded triangular" shape of 
the light curve can be an easy diagnostic for the presence 
of a comet. However, in some cases, the light curves can 
mimic an occultation by a more compact object like a clas- 
sical extra-solar planet (Lecavelier des Etangs et al. 1997). 
The color signatures of cometary occultations can be too 
small to avoid the confusion. The detection of periodicity 
appears to be critical in the diagnostic of a planetary oc- 
cultation. Alternatively, the polarization is another way to 
discriminate the two phenomena, because the light is scat- 
tered by the cometary dust roughly gathered in the same 
plane. But, with a level of at most few percents in polar- 
ization of solely the scattered light, this gives < 0.01 % 
polarization in total and would also be difficult to detect. 
In cases of planetary and cometary occultations, spectro- 
scopic follow-up observations should be planned to allow 
a better analysis of a suspected on-going detection. 



4. Probability of detection 

From the above modeling, we can evaluate the probability 
to detect a photometric stellar variation due to an extra- 
solar comet. We assume a survey of a given number of 



Fig. 3. Plot of the photometric variations due to a comet 
similar to the comet of the Fig. 2 (histogram). Here the 
result is given as a noisy observation with a photometric 
measurement every hour and an accuracy of 3a = 10~ 4 . 
For comparison, an occultation by a planet with a radius 
of 9 000 km and orbiting at 2 AU is plotted with limb 
darkening effect taken into account (thick line). The de- 
tected variation is very similar to a planetary occultation, 
the main difference being extended wings in the cometary 
case. 



stars at a given accuracy. We carried out a large number 
of simulations of cometary occultations in all directions. 

4-1. Orbital parameters 

The characteristics of the comets are given by their dis- 
tribution in the parameter space. Concerning the orbital 
parameters, the periastron density distribution is chosen 
to be the same as the one observed in the solar system: 
dn(q) cx q az dq with q G [0.1 AU, 2 AU] (A'Hearn et al., 
1995). Many sun-grazing comets have been discovered by 
SOHO with perihelion q < 0.1 AU (Kohl et al. 1997; see 
also the numerous IAU Circulars). But it is still difficult 
to infer the distribution for these small pcriastrons. The 
above distribution is thus likely biased by an underes- 
timate of the number of comets with small periastron, 
because such comets are difficult to observe. These star- 
grazing comets give a larger photometric variation. There- 
fore, we possibly slightly underestimate the probabilities 
of detection. 

We fixed all the apoastrons at 20 AU, the longitudes 
of periastron, the ascending nodes and the inclinations are 
chosen randomly. 
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Fig. 4. Same as Fig.l. Here all the parameters are the 
same except that the periastron is at 0.3 AU and the par- 
ticle size distribution is taken with a p = 0.25/xm. The 
production rate corresponds to 2 • 10 7 kg s _1 at 0.3 AU. 
Because there are more particles smaller than the wave- 
length, we see a dipper occultation in the blue by about 
15%. 



4-2. Dust production rate 

One of the most important parameter is the dust produc- 
tion rate. Because, it dictates the mass of dust in the tail, 
and constrains the detectability of the comet. The ampli- 
tude of the photometric variation is roughly proportional 
to that parameter. 

We consider that the dust production rate (P) is pro- 
portional to the area of the comet's nucleus, and that 
the distribution of the comets' radii (R c ) is similar to 
the distribution observed for the comets, asteroids, and 
Kuiper belt objects of the solar system. We assume that 
the number density of objects with radius in the range R c 
to R c + dR c is dn(R c ) cx R^ y dR c . 7 is a positive number, 
typically in the range 3 to 4 (Luu 1995). It is constrained 
to be ~ 3 - 3.5 for the comets' nuclei observed with 
R e [0.1km, 100km] (Fernandez 1982, Hughes & Daniels 
1982, Brandt et al. 1997) The same distribution is consis- 
tent with the observation of the Kuiper belt objects: 7^3 
with R e [100 km, 400 km] (Jewitt 1996), or with theoret- 
ical models for the formation of these objects: 7 ~ 4 for 
R & 200 km (Kenyon & Luu, 1998). With the assumption 
that P cx R%, the number of comets with a production 
rate between P and P + dP is dn(P) cx p-^+^^dP. We 
adopt 7 = 3.5. Changing 7 to 3 or 4 would change the 
probability of detection by less than a factor of two. 
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Fig. 5. Plot of the number of detections of cometary oc- 
cultations as a function of the photometric accuracy. Two 
cases are considered: the "pessimistic case" which is a sur- 
vey of n s tar = 30000 sun-like stars, during T a b s = 1 year 
(n CO met = 100, Pq = 10 3 kg s _1 ), and the "optimistic 
case" where among the 30000 stars, 1% (n star = 300) have 
a (3 Pictoris-like activity (n comet = 100, Po = 10 6 kg s^ 1 ). 
With an accuracy of 10~ 4 , about 10 to 10 3 detections can 
be expected. The large difference between the two cases 
shows that this kind of survey will also give information 
on the planetary evolution. 

For comparison, the dotted histogram gives the number of 
detection of planets assuming that each star has a plane- 
tary system like the solar system. With an accuracy larger 
than 10~ 4 , this gives the number of detection of giant 
planets, and mainly Jupiter-like planets. With an accu- 
racy better than 10~ 4 , Earth-like planets will be detected. 
Each step in the histogram represents the possibility to 
detect successively the Earth (8 • 10~ 5 ), Venus (7 • 10~ 5 ), 
Mars (2 • 10~ 5 ) and Mercury (10~ 5 ). We see that accu- 
rate photometric survey should detect more comets than 
planets. 



Finally, this distribution is normalized by n comet , the 
number of comets per unit of time passing through the 
periastron with a production rate larger than P . Hence, 
we have n comet cx dn(P). Thus, 



^cornet ^ Pq 



4-3. Results 



(l-7)/2 



(12) 



Using a large number of various comets, we calculate the 
probability of detection at a given photometric accuracy. 
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Then, the number of possible detections is simply this 
probability multiplied by n comet , the duration of the ob- 
servation T b s and the number of stars surveyed n sta r- 

We suppose that the time scale between each measure- 
ment is small enough 1 hour) that each variation above 
the detection limit will effectively be detected. As an ex- 
ample, we take T a b s = 1 year and n sta r = 30 000, which 
is the order of magnitude for the future space mission 
COROT. 

The expected number of detections is plotted in Fig. S. 
We considered two types of planetary systems. The first 
one is similar to the solar system with n come t = 100 comets 
per year with P > P = 10 3 kg s _1 at r = 1 AU. We see 
that few dozens of comets could be detected at an accuracy 
of 10~ 4 . We consider it as the pessimistic case. 

The second type of planetary system is supposed to be 
a young planetary system with a large cometary activity as 
during the youth of the solar system. The typical example 
is the well-known star j3 Pictoris where comets' infalls are 
commonly observed through UV and optical spectroscopy 
(see e.g., Ferlet et al. 1987, Lagrange et al. 1988, Beust 
et al. 1990, Vidal-Madjar et al. 1994, 1998). For such a 
planetary system, n comet = 100 comets per year with P > 
Pa = 10 6 kg s" 1 at r = 1 AU (Beust 1995, Beust et 
al. 1996). Note that (3 Pictoris is young but on the main 
sequence (Crifo et al. 1997), its age is few percents of the 
age of the solar system. Therefore, in a set of 30 000 stars 
there should be about n star ~ 300 stars with about the 
same activity as f3 Pictoris. Thus, few thousands comets 
could be expected with a survey at 10~ 4 accuracy. We 
consider it as the optimistic case. 

The solar system is likely not exceptional. The bottom- 
line in Fig. U is a good estimate of the lower limit of the 
expected number of detections. Younger stars may have a 
higher level of activity, with a larger number of comets; 
but, although infrared excess have been observed around 
many main sequence stars, (3 Pictoris is certainly a very 
peculiar case (Vidal-Madjar et al. 1998). Thus, the top- 
line in Fig. H gives a good estimate of the upper limit of 
the expected number of detections. 

Note that the major assumption that the dust pro- 
duction rates in the solar system can be extrapolated to 
large values is realistic. A large rate has effectively been 
observed in the recent comet Hale-Bopp where it reached 
of few times 10 5 kg s _1 at about 1 AU (Rauer H. et al., 
1997, Schleicher et al., 1997, Senay et al. 1997, Weaver et 
al., 1997). 

It is very likely that in the near future a large number 
of extra-solar comets will be detected through occulta- 
tions. 

For comparison, we also evaluate the probability to de- 
tect planets assuming that each star has a planetary sys- 
tem like the solar system. With an accuracy better than 
10~ 2 , Jupiter can be detected. Below 10~ 3 , other giant 
planets are also detectable. But because of their large or- 
bital periods, their contribution to the number of detec- 



tions is smaller than the one from Jupiter. With a total 
probability of ~ 10~ 4 to detect a giant planet in one year, 
the survey of 30000 stars will give ~ 3 planets. With an 
accuracy better than 10~ 4 , Earth-like planets start to be 
detectable. Because of their smaller distance to the star, 
they have larger contribution to the probability of planet 
detection. The Earth and Venus can be detected with a 
probability of ~ 10~ 2 in one year. Mars and Mercury, 
which are detectable with an accuracy of ~ 10" 5 , give a 
total probability of ~ 3 • 10" 2 . The comparison with the 
number of detection of comets shows that accurate pho- 
tometric surveys should detect more comets than planets. 

5. Conclusion 

We performed detailed numerical simulations of stars' oc- 
cultations by extra-solar comets. We extracted the appar- 
ent photometric variations of the central stars due to these 
putative comets. We have shown that: 

1) Extra-solar comets can be detected through pho- 
tometric variations due to occultation by dusty tails. 
In many cases, the light curve shows a very particular 
"rounded triangular" shape. However, in some remaining 
cases, the curve can mimic a planetary occultation. 

2) The photometric variations due to cometary occul- 
tations are mainly achromatic. This property will allow 
to discriminate the occultations by comets from intrinsic 
stellar variations. However, the confusion with planetary 
occultations cannot be efficiently removed by color mea- 
surements in the optical. 

3) The number of detections which can be expected 
from a large photometric survey of several tens of thou- 
sand of stars at high accuracy (10~ 4 ) is of the order of 
several hundreds of occultations per year. 

These detections will allow to explore the evolution 
of the cometary activity through the correlation with the 
stellar age. Indeed, we know from the solar system ex- 
ploration that the cometary activity significantly changed 
with time. Moreover, comets are believed to be the primi- 
tive bricks of the planetary formation, and planets per- 
turbation are needed to push them toward the central 
star. It is thus clear that the detection and analysis of 
the cometary activity around nearby stars will give im- 
portant information on the structure and evolution of the 
planetary systems. 
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